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Abstract. - A self-consistent variational theory for globules based on the uniform expansion 
method is presented. This method, first introduced by Edwards and Singh to estimate the 
size of a self-avoinding chain, is restricted to a good solvent regime, where two-body repulsion 
leads to chain swelling. We extend the variational method to a poor solvent regime where the 
balance between the two-body attractive and the three-body repulsive interactions leads to 
contraction of the chain to form a globule. By employing the Ginzburg criterion, we recover 
the correct scaling for the ^-temperature. The introduction of the three-body interaction term 
in the variational scheme recovers the correct scaling for the two important length scales in the 
globule — its overall size R, and the thermal blob size £t- Since these two length scales follow 
very different statistics — Gaussian on length scales £t, and space filling on length scale R — our 
approach extends the validity of the uniform expansion method to non-uniform contraction 
rendering it applicable to polymeric systems with attractive interactions. We present one such 
application by studying the Rayleigh instability of polyelectrolyte globules in poor solvents. 
At a critical fraction of charged monomers, f c , along the chain backbone, we observe a clear 
indication of a first-order transition from a globular state at small / to a stretched state at 
large /; in the intermediate regime the bistable equilibrium between these two states shows the 
existence of a pearl-necklace structure. 



Introduction. - The uniform expansion method of Edwards and Singh is a well known 
self-consistent variational approach to study the size and the probability distribution of the 
end-to-end distance of a chain in a good solvent [1,2]. The method is based on the uniform 
expansion of a chain in terms of the expansion in an unknown step length b\ such that the size 
of the self-avoiding chain consisting of N segments is governed by Gaussian statistics, R = 
N x / 2 bi. The variational parameter b\ is then determined self-consistently by a perturbative 
calculation, usually truncated at first order. In particular, the method recovers the correct 
scaling, (R 2 ) w b 2 (v/b 3 ) 2 / 5 N 6 / 5 , for the size of an excluded volume chain, a result which is 
consistent with the classical Flory theory [3] . In contrast to other analytical approaches for 
instance, the self-consistent field theory by Edwards [4], de Gennes [5], and the more rigorous 
renormalization group treatment by Freed [6, 7] — the advantage of the uniform expansion 
method is its mathematical simplicity in treating the complex excluded volume problem, 
which renders it applicable to more complicated systems with additional interactions. It has 
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become a standard method in polymer physics to understand the equilibrium behaviour of 
wide variety of polymeric systems. 

The uniform expansion method though applicable to wide variety of problems is only 
restricted to the behavior of polymers in a good solvent regime, where the two-body repulsive 
interaction swells the chain. As the temperature is reduced from the 9 temperature to a 
poor solvent regime, the two-body interaction becomes attractive; therefore, any attempt to 
extend this approach to a poor solvent regime requires the stabilizing influence of the three- 
body repulsive interaction. In a simple scaling form, the free energy is given by 

a «- R 2 N 2 N 3 

where v ~ rb 3 and w ~ b 6 are the strength of the two and three body interactions respectively. 
The reduced temperature r describes the distance from the temperature, i.e., r = \T— 6|/9. 
Close to the temperature, v w 0, and the chain is essentially Gaussian, R « N x / 2 b] the 
temperature can thus be determined by the balance of the first two terms in the above 
equation, which yields uiV 2 /(iV 1 / 2 6) 3 « 1 implying r « A -1 / 2 . The latter is the Ginzburg 
criterion which determines the temperature at which the entropic effects (fluctuations) become 
less important. As the solvent quality becomes poor, the polymer conformation is determined 
by the balance of the two-body attraction and the three-body repulsion; the comparison of the 
two opposing interactions (which amounts to minimizing the above free energy with respect 
to R and then comparing the second and the third term) gives rise to a globular structure of 
size R w 6r _1 / 3 A 1 / 3 , where b is the step length. Thus the size of the globule depends not 
only on the number of segments, A, but also on the reduced temperature, r. 

Apart from the overall size of the globule there is another important length scale in the 
problem — the correlation length of the density fluctuations which defines the thermal blob size 
£t [8]. For length scales smaller than £ T , the chain is unaffected by the attractive interaction 
and essentially retains the Gaussian statistics, i.e., £r w bgT 1 ^ 2 , where gx ~ 1/t 2 is the 
number of monomers in the thermal blob of size £t ~ b/r. On length scales larger than £t 
the chain begins to feel the attractive influence of the two-body interaction and the thermal 
blobs are space filling, that is, p£j, « 1; the overall size of the globule is given by R « Ct^t 1 ^ 3 , 
where nx = A/ gr is the number of thermal blobs in a globule. In contrast to the "uniform" 
expansion in a good solvent regime, the presence of the two different length scales in a globule 
require a "non-uniform" contraction of the Gaussian chain. In what follows we extend this 
variational method to a poor solvent regime by including the repulsive three-body interaction 
along with the attractive two-body interaction term. 

Let us first review some of the main results of the uniform expansion method, which is 
essentially based on defining a new step length, &i 3> b, such that the mean square end-to-end 
distance of the chain in presence of the excluded volume is Gaussian and is governed by (R 2 ) = 

Nb\. This condition amounts to replacing the original Hamiltonian Ho — ^ J Q N dsr(s) 2 with 
the reference Hamiltonian Hi = tJt J q N dsr(s) 2 . This step can easily be carried out by adding 
and subtracting Hi from Ho to give the following expression for the mean square end-to-end 

/ R 2N _ / V[r(s)]R 2 expKtti + (Hp - Hi) + H 2 )} 

^ I JV[r( s )}ex P [-(Hi + (Ho-Hi) + H 2 )} 1 U 

where Hi — J Q N ds J Q N ds'5[r(s) — r(s')] represents the two-body interaction, and v = rb 3 is 
the strength of the interaction. The idea is to expand (R 2 ) in a perturbative series about the 
reference Hamiltonian Hi such that to the first order correction in v one obtains the following 
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variational equation for the unknown parameter b\\ 

(R 2 ) 1 {Ho -Hi+n 2 ) 1 - (R 2 (w - Wi + n 2 )) 1 = 0, (3) 

where (• • • } denotes the average with respect to the probability distribution 



V[r(s)] oc exp 



26? Jo 



dsr(s) 



The Gaussian nature of the probability distribution makes the evaluation of the above averages 
simple. The details of the calculation can be found in References (1) and (2); here we simply 
present the result: 

, 6\ 1/2 N 1 ' 2 

or - or 



where a = R/Ro is the expansion factor and Rq — N 1 / 2 b. For a swollen chain, a ^> 1, the 
above equation can be solved to give R « |w| 1 /552/5 iV 3/5 ; which 

recovers the correct scaling 

for the excluded volume chain. It is interesting to note that inclusion of the higher order terms 
only changes the numerical factors without altering the overall structure of R [1,2]. 

In presence of the attractive two-body interaction term Eq. has to be modified to 
include the effects of the three-body repulsive interaction. The resulting expression is given 

by 

(r 2 ) 1 (Ho -n 1 -n 2 + n 3 ) 1 - (n 2 (n -Hx-n 2 + n 3 )) 1 = o. (5) 

where w = b e is the strength of the three-body interaction and 

H 3 = T ds ds> / ds"6[r(s) - r(s')]S[r(s') - r(s")]. 
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The quantity we are interested in evaluating is (R 2 ) x (^3)1 — (R 2 ^)^ where, as before, the 
average with respect to the Gaussian distribution function can easily be performed to give 

nih A r°° r°° r N r s r s ' 

(K 2 ) 1 {H 3 ) 1 -{ll 2 n 3 ) 1 = / dq dk ds ds' ds" [(s~s') 2 k 4 q 2 

y l /7r J Jo Jo Jo Jo Jo 

+ (s' - s ") 2 k 2 q 4 }e- k2b ^ s ~ s '^ 6 e- q2b ^ s '- s '">^. (6) 

It is clear from the structure of the above expression that the integration over k and q fol- 
lowed by integrations over the internal coordinates of the chain (e.g., s" , s' and s) leads to 
a divergence. This divergence is a common feature of excluded-volume problems in polymer 
physics, and has it origins in the use of delta function pseudopotentials. This divergence, 
though present even in the case of the two-body interaction, cancels in the overall expres- 
sion for (R 2 ) (see References (1) and (2) for more details). This is not the case when the 
three-body interaction is included, and the evaluation of above integrations turns out to be 
non-trivial. However, at the present level of the mean-field description, it is possible to make 
reasonable approximations to probe the relevant length scales in a globule. 

The idea is to first evaluate the k integration in the first term and the q integration in 
the second term of Eq. JHJ followed by integrations over the internal coordinates of the chain 
(s" , s' and s). This order of integration is important to be able to perform all but the last 
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integration without encountering any divergence. The result of such an operation can be 
expressed as 



2\ 3/2 r 1/2 7T 3 / 2 



P) 



which diverges as q — > 00. However, within the mean-field description, the divergence can 
be removed by introducing an upper cut-off for the wave number q. This is not a serious 
approximation since density fluctuations in a globule are extremely small, and the mean-field 
description suffices. Thus by introducing an upper cut-off, q m ax = 27r/A, one can probe the 
relevant length scales, A, in a globule. The result of the above integration when substituted 
into Eq. (0, produces the following equation for a: 

a 5 -a 3 +TN 1 / 2 -^^ = 0, (8) 
or A 

where we have used v = rb 3 and w = b 6 . For convenience, we have ignored all the numerical 
coefficients. First, we begin by asking for the solution of the above equation for length scale 
of the size of the chain, A = R. Since the only known length scale for the average size of the 
chain is R = N x ^ 2 bi, it is reasonable to substitute A » N x / 2 bi into the above equation to 
yield 

a 5 - a 3 + tN 1 ' 2 - 4t = 0. (9) 
or 

In presence of the attractive potential, the chain adopts a globular conformation, and entropic 
effects become less important. We can, therefore, use the Ginzburg criterion to estimate the 
temperature by comparing the first term (containing effects due to entropy) with the second 
term (containing effects due to the attractive potential) in the above equation for the condition 
a < 1. Such a comparison yields r w N^ 1 / 2 , a result which is consistent with the scaling 
analysis. The balance between the second and the third term, on the other hand, yields the 
size of the globule, R = br^^N 1 / 3 . 

Second, for length scales smaller than the size of the globule but larger than the step length 
b, we look for the density correlations of the order of the size of the thermal blob, i.e., A » £t- 
The substitution of the latter expression into Eq. J5J) yields 

/V 1 / 2 ?, 

a 5_ a 3 +rjV l/2_^L & _ (10) 

cr£r 

Since the statistics inside a thermal blob is Gaussian, we determine £t for the condition a = 1 . 
The latter condition produces £t = b/r, which, as discussed in the Introduction, represents 
the size of the thermal blob. Moreover, the Gaussian statistics for £t implies £t = g^^b, 
where qt is the number of monomers inside the thermal blob. The latter result along with 
£t = b/r leads to gr = 1/t 2 - When the size of the thermal blob £t is of the order of the step 
length b, we recover the condition, r ~ 1, defining the most compact globule of size given by 
R w bN 1 ' 3 . 

In contrast to the scaling theories, the present approach provides an alternative way to 
estimate the blob size. In presence of the repulsive three-body interaction, the resulting varia- 
tional equation needs to be regularized, a requirement that leads to a nonuniform contraction 
of the chain. The presence of the three body interaction is, therefore, essential to compute 
the overall size as well as the internal structure of the globule. 
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Thus our extension of the method of Edwards and Singh to the nonuniform case recovers all 
the essential features of a globular structure, and yields a simple variational scheme which can 
be applied to systems with attractive interactions. In the next section we use this approach 
to study the Rayleigh instability of charged globules in poor solvents. 

Poly electrolytes in poor solvents. - In addition to the two-body attractive and three-body 
repulsive interactions, the conformation of polyelectrolytes in poor solvents is also governed 
by the long-range electrostatic interaction between charged monomers. With the increase in 
the fraction of charges, /, along the chain backbone, a conformational transition occurs from a 
compact globular state at small / to a highly extended rod-like state at large /. This transition 
is due the repulsive nature of the Coloumb interaction, which at a critical fraction of charges 
makes the large globular structure highly unstable; the system gains energy by splitting into 
smaller charged globules connected by a narrow string to form what is commonly known as 
a pearl-necklace — a string of locally collapsed pearls (globules) [9, 10]. This instability of a 
charged globule is believed to be similar to the Rayleigh instability of a charged liquid droplet 
which occurs as soon as its electrostatic energy, Q 2 k B Tl B /e 2 R, exceeds the surface energy, 
7-R 2 , i.e., Q > (7i? 3 /fcsTZs) 1 / 2 e [11, 12], where Q is the charge on the droplet of size R, e 
is the elementary charge, Ib — e 2 / A-neksT is the Bjerrum length, e is the dielectric constant 
of the solvent and 7 is the surface tension. As a result of this instability a charged droplet 
spontaneously breaks up into smaller droplets, placed infinitely apart, each with a charge 
less than the critical one. However, due to the constraint of chain connectivity the preferred 
equilibrium conformation of a charged globule is that of a pearl-necklace, which with the 
increase in / results in a stretched state. Although scaling theories do provide some evidence 
of the abrupt transition at a critical fraction of the charge from a globule to a stretched state, 
a clear indication of its being first-order is still missing. In what follows we study in detail 
the nature of the transition that brings about the Rayleigh instability in a charged globule; 
the approach we use is the variational method discussed in the last section. 

The repulsion between charged monomer is described by the Coulomb interaction given 
by Tic = Tf f ds L ds'\r(s) — r(s')| -1 , where v c = J 2 Ib is the strength of the Coloumb 
interaction. This electrostatic interaction when included into Eqs. J2J) and (JSJ yields the 
following equation: 

(R 2 ) 1 {Ho -H 1 -H 2 + H 3 + -He), - (R 2 (W -H 1 -H 2 + H 3 + H c )) 1 = 0. (11) 

The additional term that needs to be calculated is given by (R. 2 ^ (W c )i — (R 2 W C ) 1 . The 
details of the calculation will be presented elsewhere [14]; here we simply present the result: 

(R 2 ) 1 (H c ) x {&H c ) l = -± ^) 1/2 fl B N^h (12) 

The above expression when substituted into Eq. pi|l gives the following equation for a: 

a 5 -a 3 + tN 1 ' 2 - -L - ( N 3 / 2 a 2 = 0. (13) 
a 6 \ J 

The above equation for a can be solved in the following two limiting cases: when f 2 lBN/b <C t, 
the attractive influence of the two-body interaction dominates, and the comparison of the 
third and the fourth term yields a globule of size R = br^ 1 / 3, N 1 ^; in the opposite limit of 
plsN/b ^> r, the repulsive effects of the electrostatic interaction overcomes the attractive 
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Fig. 1 - The numerical solution of Eq. I I -'it for a is plotted as a function of the strength of the 
electrostatic energy / = f 2 lsN/b for r = 0.5 and N = 1000. The region between the dotted lines is 
the region of coexistence between a globule and an extended state. The bistable equilibrium between 
these two states is represented by the dashed line, where R indicates the rod-like extended state 
and G indicates the globular state; the unstable part of the region is represented by the curve U. 
f c represents the critical fraction of charge at which the first-order transition from a globular to an 
extended state occurs. The parallel horizontal lines represent regions Ai and A2 of equal area. 



interaction, and the comparison of the first and the last term yields a — (/ 2 ^s/6) 1 / 3 7V 1 / 2 im- 
plying a stretched state, R — 6iV(/ 2 /s/6) 1 / 3 . In terms of the electrostatic blob picture the lat- 
ter represents the linear arrangement of N/g e i (electrostatic) blobs of size £ e ; = (/ 2 Zb/&) -1 / 3 , 
where g e i — (jBf 2 /b)~ 2 ^ 3 is the number of monomers in a blob. Therefore, the expressions 
obtained from the variational approach are in complete agreement with the scaling predic- 
tions [9,10]. 

In the intermediate regime, / = f 2 lsN/b ~ r, Eq. i|13fl can only be solved numerically. 
The result is presented in the form of Fig. (1), which is a plot of a as a function of the strength 
of the electrostatic energy. In the intermediate regime (the region between the dotted lines 
represented by f m in an d fmax) the curve folds back on itself to give a bistable equilibrium 
between a globular state (a < 1) and a rod-like state (a > 1). The coexistence of these two 
states in the region between the two dotted lines in Fig. 1 is a clear indication of the existence of 
the pearl-necklace structure in this region. Moreover, the abrupt transition from a globular to 
an extended state is evidently first-order in nature. Therefore, we use the Maxwell equal-area 
construction to determine the critical fraction of charge f c at which the spontaneous Raylcigh 
splitting of the charged globule occurs. The point f c in Fig. 1 is the point which ensures that 
the area under the two curves (represented by A\ and A2) is equal. The intermediate regime 
corresponds to the balance between the forces close to the Rayleigh splitting, i.e., when the 
electrostatic force /csT/s/ 2 iV 2 /i? 2 n becomes equal to the line tension ksTr/b. The chain 
extension then scales as R pn = {IbP /6r) 1 / 2 iV6, where R pn represents the mean size of the 
pearl-necklace structure. 
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The pearl-necklace structure, as seen in simulations, is a highly fluctuating structure both 
in terms of the position and the size of pearls [13]. In the mean-field description, the existence 
of such a structure is evident through a bistable equilibrium between a globular and a stretched 
state, as is clear from the coexistence region, f m %n < / < fmax, in Fig. 1. The scaling theories 
usually estimates the average size of a pearl-necklace by comparing its surface energy with 
that of the electrostatic energy. In a future publication, we will present a detailed analysis 
of the dependence of R, in the coexistence region, on /, N and r to make comparison with 
the scaling predictions. The effects of the variation of r and N on the coexistence region 
(pertaining to pearl-necklaces) and f c will also be discussed [14]. 

Conclusions. - We have extended the uniform expansion method of Edwards and Singh 
to a poor solvent regime by accounting for the three-body repulsive interaction in addition 
to the two-body attractive potential. The self-consistent variational approach recovers the 
two important length scales in a globule — its overall size and the thermal blob size. This 
method can be applied to several different polymeric systems where the effective two-body 
potential is negative. To illustrate the generality of our approach, we applied this method to 
polyelectrolytes in poor solvents to study the Rayleigh instability of charged globules in detail. 
The presence of the additional long-range repulsion due to the electrostatic interaction between 
the charged monomers destabilizes the globular structure as the fraction of charged monomers 
along the chain backbone is increased. At a certain critical fraction of charges the globule splits 
into locally collapsed monomers joined by a linear string to form a pearl-necklace. Using the 
present approach, we observe a clear indication of a first-order transition, at a critical fraction 
of charged monomers, from a globular to an extended state. The intermediate regime of 
bistable equilibrium between these two states is the region where the pearl-necklace structure 
is stable. 

In contrast to flexible polymers, the behaviour of semiflcxiblc polymers in poor solvents is 
not well understood. Given the simplicity of the present approach, it holds out the possibility 
of addressing many such issues, including the equilibrium behaviour of biopolymers with 
specific interactions [14]. 
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